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Abstract 

A recently proposed set of sum rules, based on the pion-Kaon scattering amplitudes 
and their crossing-symmetric conjugates are analysed in detail. A key role is played by 
the I = 7T7T — > KK amplitude which requires an extrapolation to be performed. It is 
shown how this is tightly constrained from analyticity, chiral counting and the available 
experimental data, and its stability is tested. A re-evaluation of the 0(p ) chiral couplings 
Li, L2, L3 is obtained, as well as a new evaluation of the large N c suppressed coupling L4. 
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1. Introduction 



Chiral perturbation theory is a rigorous approach to QCD in a restricted but nonper- 
turbative regime, which has recently been developed to 0(p 6 ), i.e. to the next-to-next to 
the leading order Foundations of this method Q and the abundant work which has 
followed the basic papers where the NLO theory was defined [||, || are summarized in the 
review [||. The chiral expansion is based on an effective field theory and, as such, involves 
an increasing number of coupling constants with increasing chiral order. In SU{3), ten 
couplings Li(n) are involved at 0(p ) and ninety more couplings Ci(fi) appear at the next 
order. In order to make predictions to 0(p 6 ) accuracy, estimates of the Ci(fi) must be 
performed but, moreover, the values of the Li may have to be modified, compared to their 
determination using 0(p 4 ) accuracy. The order of magnitude of such a variation, which re- 
flects the rate of convergence of the expansion in the strange quark mass, can be estimated 
by comparing several different 0(p 4 ) determinations of the same coupling constants. This 
is one purpose of the present work in which we propose a new determination of Li, L2 and 
L3 from a set of sum rules based on the pion-Kaon amplitude and its expression in ChPT 
at one-loop ||. We will compare these results with the previous determination from the 
K14 form- factors |7|, || and the (partial) determination from tttt sum rules || . 

A few 0(p 4 ) coupling constants are still very poorly known, in particular, those which 
are suppressed in the large N c limit: L4 and Lq. Naively, one may even question whether 
such a suppression should actually hold, because these couplings were shown to be con- 
trolled by physics of the scalar meson resonances § which fail to obey simple large N c 
rules. On a more sophisticated level, one may note that some of the large N c suppressed 
mechanisms, like internal quark loops, are partly taken into account in the chiral expansion 
via meson loops. The question remains of what value of the scale \i is the one at which 
the suppression operates. Another related interesting issue is that of the phase structure 
of QCD-like theories as a function of the number Np of massless flavours and the value of 
Np lt above which chiral symmetry is no longer spontaneously broken. Some recent lattice 
simulations [1C] have obtained values as small as Np 11 ~ 4. If true, this should affect the 
SU(3) chiral expansion. For instance, it can be seen that L4 and Lq control how the chiral 
order parameters F n and < uu > respectively evolve from Np = 2 to Np = 30. Clearly, 
a small value of N F rit , should lead to anomalously large values of L4, L§. 

In view of this, an interesting outcome of the present work is a determination of L4. In 
principle, it could have been extracted from the K14 form-factors, but this is not feasible 
in practice because its contribution is accidentally suppressed [0, [|. Here, we will take 
advantage of the fact that no such suppression affects the ttK amplitude and we will show 
that an evaluation of L4 is then possible, which is at the same level of reliability and 
accuracy as that of L±, L2, L3. Several recent papers have considered aspects of the 
pion-Kaon scattering amplitudes [11, 12, |DJ. One purpose is a better understanding of 
the scalar resonances (see e.g. []l4| for a recent discussion of the experimental situation). 
This question, of course, is not unrelated to that of the size of the chiral couplings M. 

The dispersive formalism on which the sum rules are based has been developed in 



a previous paper [15]. This formalism is reviewed in sec. 2 below and presented in a 
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form suitable for comparison with the 0(p 4 ) expression of the amplitude, which has been 
computed some time ago by Bernard et al. ||, as well as the 0{p & ) expression which 
should be available in the near future. The detailed form of the sum rules for the 0(p 4 ) 
coupling constants are then presented in sec. 3. The practical evaluation of these sum 
rules, making use of the available experimental data then presents a difficulty: because of 
s — t crossing the irir — > KK amplitude appears and it is needed below the experimentally 



accessible energy range. This was noted in ref. [15| in which only results not depending on 
this amplitude were presented. Extrapolation of the irir — > KK amplitude, in particular 
for the S-wave, is a problem which was considered a long time ago [jl6| , |l7fl . We discuss 
this question in some detail in sec. 4, and then present all the results. 

2. Dispersive representation, crossing-symmetry and chiral counting 

Basic work on dispersion relations related to the pion-Kaon amplitudes has been re- 
viewed by Lang pq|. In order to determine the number of subtractions we make the 



assumption that standard Regge phenomenology applies. As shown in ref. [15] the disper- 
sive representation can be recast in a specific form by taking into account chiral counting. 
Dropping terms which are of chiral order 0(p 8 ) it can be put in a form which involves 
functions of only one of the Mandelstam variables s, t, u and are analytic, except for 
a right-hand cut, plus a polynomial. This was first demonstrated for the case of the 
pion-pion amplitude in ref. |19[ |. Let us begin by recalling some basic facts and some 
notation. 

2.1 Notation and conventions 

Making use of s — u crossing, the two independent isospin / = 1/2 and 1 = 3/2 
pion-Kaon amplitudes can be expressed in terms of the 1 = 3/2 one, 

1 3 3 3 

F^(s,t,u) = -~Fl(s,t,u) + -Fl(u,t,s) . (1) 

It is convenient to introduce the amplitudes F + and F~ which are respectively even and 
odd under s — u crossing because they require a different number of subtractions. In terms 
of isospin amplitudes, they are defined as 

F+(s,t,u) = -F5(s,t,u) + -F5(s,t,u) 

F~(s,t,u) = ^F3{s,t,u) - ^F-2(s,t,u) . (2) 

Under s — t crossing, one generates the 1 = and 1 = 1 irir — » KK amplitudes, 

G°(t,s,u) = V6F + {s,t,u) 

G l (t,s,u) = 2F~(s,t,u) . (3) 
The partial wave expansion of the irK isospin amplitudes are defined as 

F / ( S ,t) = 167r^(2/ + l)P i (z s )//( S ) . (4) 
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In a similar way we can expand F + and F , the corresponding partial-wave projections 
are denoted f^~{s) and fr(s). The s-channel scattering angle appearing above is given by 

s(t — u) + m 2 _m\ 



with m± = rriK ± "lyr • 



(5) 



s (s — m 2 _)(s — m+) 
The partial- wave expansion of the tttt — > i^i^ amplitude is conventionally defined as, 

G 7 (t, a) = 16vrV2 ]T(2Z + l)P^)<?f (t)(q„(t)q K {t)) 1 . (6) 



with 



Ami 



lp > Zt ~ Aq^t)q K (t) ■ [<) 

The rationale for introducing the factor {q n qK) 1 in eq. (§) is explained by Frazer and Fulco 
[poll- It ensures that the partial- wave amplitudes <7jf(i) have good analytic properties. With 
these definitions, the partial-wave S-matrices are given by 



J (s - m 2 _)(s - mV) 
ttK^ttK: Sf(s) = 1 + 2i±l —fi(s) 



tttt -» KK : Sf(t) = Ai 



2.2 Dispersive representation of F + {s,i) 



Mt)qK(t)) l+1/2 i 



(8) 



One first writes down a dispersive representation with t fixed (and small). According 
to Regge phenomenology, the asymptotic dependence as a function of s is controlled by the 
Pomeron, implying the need for two subtractions, which would also result on the general 
basis of the Froissart bound, 



r,+ , x /x I f°° ds' f s 2 



+ 



ImF + (s',t), 



(9) 



giving F + (s,t) in terms of an unknown function of t. Next, following ref. [19], one splits 
the integration range into two regions a) [m\, A 2 ] and b)[A 2 , oo], A being the scale of the 
chiral expansion, i.e. A ~ 1 GeV. In the lower integration range, we can apply the chiral 
counting and drop the imaginary parts of the partial waves with I > 2 which are 0(p s ), 
i.e. we put 

s'(t-u') + A 2 K7T 



ImF + (s,t) = 16tt 



/m/ + + 3/m/ 1 + (s^ 



(s' — m 2 _){s' — m\) 



s'<A 2 . 



(10) 



In the region b) we can expand in terms of s, t, u divided by A 2 again dropping terms 
which are 0(p s ). After some reshuffling of part a) and absorbing functions of t into c(t), 
one obtains the fixed t dispersive representation in the form 



F+(s,t) = c(t) + [w +(s) + (t - u)Wf(s) + (s~u) 



-(2us)- 



ds' 

7T J A 2 (s') 3 



1 + 3 



^ttK — 1/2 



ImF + {s',t) + 0(p 8 ) 



(11) 
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In the equations above we have introduced the notation 

^PQ = m 2 P + mQ, A P Q = mp-mQ. (12) 

The functions Wq (s), W^(s) are analytic except for a right-hand cut and are given in 
terms of the S and P waves of the pion-Kaon amplitude, 



W+(s) = 16 / -p- Imf+(s') + A 2 

i . . f A ' 2 ds' 2,1m f Us 1 ) , , 

<•> - 16s U —s v^w^r) ■ <13> 

In order to further constrain the function c(t) appearing in eq. ([IT]) we must write down 
for F + (s, t) a dispersion relation involving the cut in the t variable. A possibility is to use 
a dispersion relation with s fixed. Alternatively, one can use one with us fixed, us = b, 
(hyperbolic dispersion relation) which was shown to have better convergence properties 
[pi]]. In this case, the variables s and u are functions of t denoted s b and u b , 



t / ^ 2 



s b (t) = + y ^ K --) -b, u b (t) = b/s b {t) . (14) 

The function F + (s b , t) is an analytic function of t with a) a right-hand cut 4m 2 < t < oo, 
and b) a left-hand cut — oo < t < m 2 __ — b/m 2 ^ . In the following, we will adopt a specific 
value for 6, 

b = A 2 Kv = m 2 _m 2 + , (15) 

which corresponds to backward scattering, z s = — 1. In that case, the upper limit of the 
left-hand cut is t = 0. In the asymptotic regions t — ► dboo, it is simple to verify that the 
dominant divergence is controlled by the K* or K\ Regge trajectories, and it is therefore 
plausible that a single subtraction is sufficient in this case. The following representation 
is then obtained 17, 15], 



VT Jm\ S' \S' - Sb s' - u b 
t f 00 dt'ImG°(t',s' b ) 



where s' b = s b (t') (see eq. (p|) )and 



(16) 



4 = - s' - . (17) 

s' 

Next, one splits the integration range as before and drops terms which are 0(p 8 ). Equating 
the representations ([ll]) and (|l6]) then determines the unknown function in the former 
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expression leaving just one undetermined constant. Introducing the following notation for 
the various high-energy integrals which are involved, 



H + (n) 
H + (n) 



1 



ds' 



7T J A 2 {s'Y' 

ds 1 



1 



7T J A 2 (s') n 

i r°° ds' 



ImF + (s',0) 



d t ImF + (s',0) 



n 



Gt(n) 



7T J A 2 («') 

i r°° dt' 



ImF + (s',t' b ) 
ImG°(t',s' b ) , 



VQttJa 2 it') 

we finally obtain the following dispersive representation for the amplitude F + (s,t) : 

F + (s,t) = C+ [w+(s) + (t-u)W+(s) + (s~u)\ +U (t) 
A2 3Imf+(s') 



(18) 



+16* 
+2bt 



.2 S (s' — m 2 _){s' 



2us 



H + (3) + tH + (3) + 3(E nK - t/2)H + {A) 



H+{3) - 3/2F+(4)J + tG+(2) + t 2 G+(3) + i 3 G+(4) 
-i# 6 + (2) + t(i - 4£^)F+(3) - t(i 2 - Qt^ K + 12£ 2 X - 3&)i7 6 +(4) + 0(p 8 ) . (19) 

Apart from a polynomial, this expression involves the functions Wn~(z), W^(z) which are 
defined in terms of the S and P waves of the ttK amplitude (see eqs. (13)p3])) and the 
function Uq(z) which is defined in terms of the S wave of the tttt — ► KK amplitude, 



TJ < \ 16 



A2 dt' Im &w 



(20) 



This derivation shows that the specific form of the amplitude eq. (|l~9| ) must hold in chiral 
perturbation theory at 0(p 4 ) (which we will check explicitly below) and also at 0(p 6 ). 

2.3 Dispersive representation of F~(s,i) 

We proceed in the same way as for F + (s, t) by first writing down a dispersion relation 
with t fixed, the only difference is that now, the behaviour at large s is dominated by the K* 
and K2 Regge exchanges and, therefore, a dispersion representation with no subtraction 
should converge, 



F-(s,t) = - I™ ds' (V ±-)l m F-(s',t). 

7T J ml \S' - S S' — Uj 



(21) 



As before, one splits the integration range into two pieces and in the lower energy range 
one retains only the S and P waves, obtaining 



F-(s,t) = W (s) - W (u) + (t- u)W{{s) - (t - s)Wr(u) 



+ (s-u) 



16 



A 2 



3/m/f (s')ds' 



+ 



1 f°° ImF-(s',t)ds' 



,2 (V — m_){s' — m+) tt J a 2 (s' — s)(s' — u) 



(22) 
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The functions W (s) and W 1 (s) are exactly analogous to their + counterparts defined 
above, 

A2 ds ' (r^ f -f^ , A 2 3/m/f( S ') \ 



X7r (s'-m 2 .)( s '-m 2 h ) 



W -(s) = 16[ -*_(/ ro /-( a ') + A 

^> = '<^ (^V < 23) 

This representation has no undefined functions but convergence is ensured only for negative 
values of i. One can extend the range of validity in t, and also display the cut structure 
by combining with a hyperbolic dispersion relation. One writes a dispersion relation at 
fixed us = b for the function 

F~(s,t)/(s -u) (24) 
which is even in s — u and thus free of kinematical singularities, and one obtains 

'GHt', s ' b : 



F-(s b ,t) 1 [°° dt' 

-Im 



s b - u b 2tt J 4m 2 t'-t 



s 'b ~ < 



1 r°° dgf ImF-(s',t' b ) 
vr Jml (a' - s b )(s' - u b ) 



In the low energy region of the right-hand cut, only the P wave of the hit — > i<Ci<C amplitude 
will contribute, which generates the function 

C/ 1 (z) = 6v / 2 dt / Jmgll(t0 . (26) 

Equating the fixed t and fixed us representations gives the following equation, valid for 
t < 0, 

3Imf{(s') if 00 ,, ImF-(s',t) 



32 T V, 3 ™ M + if & 
V is' - mi Is' - mi ) 7r Ja 2 



,2 (s' — m_)(s' — m\) ir J a? (s' — s b )(s' — u b ) 

UAt) + -1 r d t' ImGl V><) + I r ds > IrnF-(s',t' b ) (2?) 
^ W + 2vr A* * (tf - t)(s' b -u' b ) + n J A 2 (s> - s b )(s> - u b ) [27) 

which relates the P- waves in the irK and the irir — > if if channels. Finally, introducing 
the following notation for the high-energy integrals, 

1 r°° r/s' 

G ' ( " ) = ^i ( t -)"-K-<) J "' Gl( '' 4K (28> 

we obtain the dispersive representation for F~(s,t), valid up to 0(p 8 ) contributions, in 
the form, 

F~(s,t) = W (s) - Wq(u) + (t - u)Wi(s) - (t - s)W{{u) + (s - up^t) 
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+ 



ds '- 



3/m/f( S ') 



+ G 6 -(2) + tG 6 -(3)+^G 6 -(4) 



{s _ „){ _ l 6 ^ _ {s ,_ m 2_ ){s/ _ m 2 +) 

+H b (2) + (2E wK - t)#r(3) + [(2S w a- - t) 2 - k]# fe ~(4) + (6 - (4)} 
+0(p 8 ) . (29) 



On the right-hand sides of eqs. (|T^)(^) the dependence on the cutoff A must cancel: we 
have verified that it does, up to 0(p 8 ) terms. The dependence upon the parameter b must 
also cancel. This gives rise to constraints among the nK and tttt — ► KK amplitudes and 
their derivatives which we have not explored. 

3. Chiral representation and sum rules 

3.1 Chiral representation at 0(p 4 ) 

First, let us recall, that at the leading chiral order, 0(p 2 ), one has 



or 



F$(s,t) = -^(4s + 3t-4S 7r ^), F§(a,t) = — f (-2s + 2^ K ) 



4 J 2 4 / 2 



The corresponding irK partial waves, are, first for I = 



/<?(«) 



1 



128tt/| 



5s - - 



/<?(*) 



64tt/^ 



-2s + 2T, 7T k) 



then for Z = 1 



fi( s ) 



1287T/2 



(s-2S^ + ^), / 1 f (s) = 0, 



(30) 
(31) 

(32) 
(33) 



while the partial waves for I > 2 vanish at this order. In the tttt — * KK channel, the I = 
and I = 1 partial waves are 



487T/2 



(34) 



At the NLO order, according to the discussion above, the ttK amplitudes must have the 
following form, 



Wq(s) + (t- u)wt(s) + (s u)l + U (t) 



F + (s,t) 
+A|t 2 + \%(s - u) 2 +/3 + t + a 



(35) 



and 



F-(s,t)= W (s) + (t-u)W 1 (s)-(s^u) +(s-u)U 1 (t) 
+(s-u)(A^ + /T) . (36) 
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Indeed, the calculation was performed in ChPT at 0(p 4 ) in ref. ||, and it is not difficult 
to recast their result in the above form. We display the explicit expressions below, which 
will be used in the derivation of the sum rules. The W functions receive contributions 
from ttK and rjK intermediate states, 



W, 



(37) 



For the Wq P q functions, one obtains from || 



0,itK 



9 9 9 2Ai XL a- 4 A 4 



19s z - 28sE wX + 12X^ - 9A^ W + 
4A 4 



+ 



K 71 



0,r)K 



3s 2 - 4sX 7rA + ^X 2 A . + Af^ + 6A A7r A Aj; - Kn (AKn^Kri 



2Ak v ^ttk) + 



4A 2 A 2 

Kit Ktj 



4A 2 A 2 

f / s Kir Kr\ j, 

^(s) 0) 



(38) 



In these expressions, Jpq(s) is the standard one-loop function j| which has the following 
dispersive representation 



Jpq(s) 



16vr 2 J(m P +m Q y- 



ds' 



V) 2 (s'-s) 



with 



^pq(s') = ( s ' - ( m P + m Q) X s ' - ( m P - m Q) ) 



Now the Wq P q functions are 



W 0, V K = W 0, V K 



(39) 
(40) 

(41) 



The last equality holds because the r\K state has isospin I = 1/2. Also Imf^ vanishes at 
0(p 4 ) (and also at 0(p 6 )) and consequently, 



Wt(s) = W l (s) . 



(42) 



The expression for the W^~ P q components is the same for PQ = irK or r/K and is given 

by 



A 2 

WZ PQ = (s- 2X P q + -^) J PQ (s) - 4A PQ J' PQ (0) . 



i,PQ — <"-'i'<} - '-»/'(.>•',"(.)>«" • (43) 

Finally, the XJ\ functions at 0(p 4 ) are 

16/ 4 U (t) = 2t(2t - ml)J n7T (t) + 3t 2 J KK (t) + 2ml(t - \m 2 K )3 m {t) , 
48/ 4 U 1 (t) = 2(t - 4 m 2)4({) + (t _ 4m^) J^(t) . (44) 
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The imaginary parts of the W and U functions at 0(p 4 ) can be recovered from the defini- 
tions eqs. (li~3|) in terms of the imaginary parts of fi (s), <?o(s) and g\ (s) and using uni- 
tarity to relate the latter to the / = 0, 1 amplitudes ttK — » ttK, r\K and tttt — > tttt, KK, r/rj 
computed at 0(p 2 ). For instance, unitarity gives 



\\fh*)\ 2 +l\fh')\ 2 +-^\fh')\ i 

6 6 Kk{s) 



(45) 



and using the 0(p 2 ) expressions 
imaginary part as in eq. (|38|). 



|33|) for the partial waves one recovers the same 



The separation in eqs. (|35|)(36) into a polynomial part and a part with cut-analytic 
functions is arbitrary: we have only required that each piece be scale independent (a 
different choice was made in ref. |15|| ) and finite. The coefficients of the polynomials are 
simple linear functions of the coupling-constants Lj(/f). Using the following notation, 



L P = log~Y, R P q 



m. 



m. 



m. 



log-f 



VTLf 



(46) 



\r- rap 

and the result of ref. B one finds for the coefficients entering the F~ amplitude at 0(p 4 ) 



1 2m 2 r 



1 



"L 3 + 



1 



L5 ~ 5127r 2 (® Lk + 5R ^K + RqK 

4 



512vr 2 



nr 



log — f + R vK + R 



m 



K 



(47) 



The coefficients entering the F + amplitude, then, have the following expression in terms 
of the Li's, 



a 



8m 2 m K 



J 7T 



{4Li + L 3 - 4L 4 - L 5 + 4L 6 + 2Lf. 



5127T 2 



1 2 

Lk - RttK + -^RKr) — g 



P + = (3~ + 



8(m 2 + m 2 K ) 



+ 



???; 



128tt 2 /4 



log 



f 4 
j - 

2 

m 2 



-2Li 



2 L3 + L4 + 512^ 



Lr + RttK + 



/ 4 A+ = 8L1 + 2L 2 + ^L 3 + ^2 [-8L W - 10L^ - 4R nK - 15] 



2L 2 + -L 3 H — ^ 

2 d 512vr 2 



-6Lk — SR-kK — Rt)K + 2 



(48) 



This completes the rewriting of the chiral formulas of ref. § in a form which allows easy 
matching with the dispersive representations. This matching generates a number of sum 
rules. For F~ , the dispersive formula has no subtraction constant, which implies that 
the two coefficients (5~ and Aj~ can be expressed as sum rules. For F + , one subtraction 
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constant remains and this implies that the three coefficients (3 + , A^ , A^" are expressible as 
sum rules while the fourth one, a + , remains undetermined in this approach. Using eqs. 
(E7|)(48) it is then easy to generate sum rule expressions for the Lj coupling constants, 
which are given in terms of (3^, Xf as simple linear combinations. For instance, L±, L2 
are given by 

f 4 1/4133 23\ 

L x = —(Xt - At + 2A7) = — L w - -L K + -R nK + -R„ K 

g v 1 2 -r 1 ; 5127r 2 ^ 3 6 8 8 12/ 

L 2 = &(2X+ + Ar) - -^-5 (-±1* - ^Lk - ^* - \r v k + ^) • (49) 

while L3 is immediately given in terms of A^~. The coupling L4, finally, is obtained from 
the following combination 

If 5 1 7 777^ 777 ^ 1 

2^ + -itW + - - + . _ * - log -f . (50) 



512tt 2 ^ " ' 4 ™ ' 4 ^ 2 ' 2(m2+m2 c ) i " & m2 / 

We observe that while the coupling L5 is present in the expression for /3~ , it appears 
multiplied by m 2 (not m 2 K ) and thus makes a very small correction to the leading 0(p 2 ) 
contribution. 

3.2 Sum rules 

The dispersive representation of the irK amplitudes eqs. ( |i~9| ) (29) contains one arbi- 
trary parameter, while the polynomial part of the chiral representation eqs. (|35|) (j36|) at 
0(p 4 ), contains six coefficients: comparing the two representations should yield five sum 
rules for these coefficients which will translate, in principle, using expressions (|48])(^7|) 
into sum rules for the five coupling constants ij(ju), i = 1,5. The explicit form of the 
sum rules are obtained by noting that differences like Wq (s) — W (s), in which W is 
computed to 0(p 4 ) accuracy, are analytic up to 0(p 6 ) contributions, 

Im(W+(s)-Wt(s)) = 0(p 6 ) . (51) 

Therefore, up to 0(p e ) terms, we can expand these differences as polynomials, 

W±{z) -Wf(z) = Af + Bfz + C±z 2 

Ui{z)-Ui{z) = u l + v l z + wiz 2 , (52) 
for I = 0, 1. We also introduce 

1 " b J m l {s'-m\){s'-ml) ■ [b6) 

The coefficients A^ 1 , £?, etc... are given as inverse moments of the imaginary parts of the 
ttK and nir — ► KK S and P waves, integrated between the threshold and A 2 , with the 



11 



chiral part being subtracted. Together with the integrals (|lq)(p8|) over the high-energy 
region, [A 2 , oo], they form the building blocks of the sum rules. 



Equating the chiral and the dispersive expressions, taking into account eqs. (|52j), we 
finally obtain the following sum rule formulas for the polynomial coefficients in the 0(p 4 ) 
chiral representation (|35|)([3C 



pr = -ir + A{ + Bo +m + 2E^C - + G 6 -(2) + H^{2) + 2^ K H^(3) 
/?+ = !+ + 3,4+ - B+ + v + 2^ K (2B+ -C+) + G+{2) 
-H+(2) + 2£ 7rK (tf+(3) - 2tf+(3)) 

A + = -|fif + \<% + w + G+(3) - ^+(3) + fl+(3) 

A? = ^ + ~C + + ~# + (3) 

= - C - + v x + G 6 - (3) - fT 6 " (3) . (54) 

The derivation and the structure of these sum rules are very similar to those which were 



proposed for tttt scattering in ref. [22|. We now discuss the practical evaluation of these 
formulas. 

4. Evaluation of the sum rules 
4.1 nK amplitudes 

We will make use of the two most recent high-statistics Kp production experiments, 



both performed at SLAC, which have determined irK amplitudes. Estabrooks et al. [23| 
have considered several charge combinations enabling them to determine separately the 
/ = 1/2 and the I = 3/2 combinations. For the isospin / = 3/2 it was observed that the 
P and D waves remain very small below y/s = 2 GeV: in our calculations we will only 
include the S-wave. A few years later the K~tt + — * K~tt + amplitude was remeasured in a 



slightly larger energy range by Aston et al. [24|. For the I = 1/2 S and P waves, we have 



performed fits of the data of Aston et al. with parametrisations in terms of Breit-Wigner 
plus background similar to those used in this reference using the 1 = 3/2 S-wave from ref. 
P3 l- hi these fits, we have imposed that the scattering lengths be equal to their values in 
ChPT ||. Relaxing this constraint, however, makes very little change in the results. For 
the partial waves 1 = 2 — 5, we have used exactly the same parametrisations as provided 
in ref. [24]]. Above y/s = 1.5 GeV, ambiguities arise in the determination of the S and 
P waves. Estabrooks et al find four different solutions and Aston et al. two. It has been 



pointed out in ref. [13] that one of these violate the unitarity bound; therefore, we have 
used the remaining one. In our sum rules, we note that the contribution from the S and 
P waves in this energy region becomes negligibly small anyway. Above the energy range 
covered by these experiments we use Regge parametrisations which we will discuss in more 
detail below. 

4.2 7T7T — > KK amplitudes 
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Let us first discuss the S-wave amplitude, 



9o(t) = \g8(t)\exp(i<j> s (t)) 



(55) 



which is a crucial ingredient in the sum rules and is needed for t > 4m^ while it is measured 
only in the range t > Am 2 K . Analyticity, as is well known 17, 25, pqj , is the key to 



performing this extrapolation. To start with, the phase of the amplitude, 4>s{t), may be 
considered as known in the whole energy region of interest. Firstly, in the region where 
Tin scattering is elastic, 4>s is identical to the tttt phase shift (modulo 7r) from Watson's 
theorem. It is now well established that, to a very good approximation, the domain where 
7T7T scattering is effectively elastic extends up to the KK threshold (see e.g. [^]). Above 
this point, <f>s (t) has been measured in experiments, we will use the two most recent ones: 
Cohen et al. 



| (who considered K + K production) and Etkin et al. [29 



|] (who 

considered K^Kg). These data are shown in Fig. 1 together with the curves which will be 
used in the calculations. One observes that the two data sets are compatible except very 



400 




Figure 1: Phase of the tttt — > KK I = amplitude. The curves are as described in the 
text. The three curves from bottom to top correspond to three increasing input values for 
the phase at the KK threshold, <j>s(t ) = 150°, 175°, 200°. 

close to the KK threshold. (It must be recalled here that experiments actually measure 
S — D interference and thus determine only the difference 4>s ~ 4>D- in treating the data 
of Etkin et al. we have used the D-wave model of ref. f3l[| rather than that used in the 
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original paper J29[|.) In the energy region y/t > 2m x, we fit the combined set of data with 
piecewise polynomials. In performing the fit, we have excluded the small energy region 
where the two data sets are incompatible and we have instead fixed the value of the phase 
at threshold 4>s(Am 2 K ), which we have allowed to vary between 150° and 220°. This may 
seem like a wide range and one could think of making use of the equality between <ps 
and the irn phase-shift at the threshold to improve on that. However, the available irir 
data points closest to the KK threshold have large error bars. In order to get exactly 
at the threshold, one needs to perform a fit and the result is uncertain because the irn 
phase-shift varies extremely rapidly in this region [^7|. One ends up with the same range 
of values as we have chosen. In the energy range yft < 0.8 GeV, the curve in Fig. 1 is the 
result from the recent Roy equations analysis of ref. |3^ j, with a{] = 0.22, o,q = —0.0444. 
Because of the new data p3| there is now a rather small uncertainty on the curve in 



this energy region [34]. In the energy region between 0.8 GeV and the KK threshold we 



use the simple interpolation formula, 

M*) = « + ( 56 ) 

in which the three parameters are determined from imposing continuity at both ends and 
continuity of the first derivative at the lower end. 

Once the phase is known, determining the modulus in the region [4m^,4m^] is a 
standard Muskhelishvili-Omnes |3^, ^] problem because satisfies the following integral 
equation 

vr J Ami "(.* ~~ t ) 

where A(t) has only a left-hand cut. A(t) can be expressed explicitly in terms of nK 
partial wave amplitudes and tttt — > KK partial waves with I > 2, 

/*oo /*oo 

A(t)=J2 ds'K m {s\t)Imf+{ S ')+ / dt'G l0 (t',t)(q n q K ) l Imgf(t') . (58) 

I Jm % l>2,even J4m * 

For t < 1 GeV 2 , A(t) is dominated by the ttK S and P waves and eq. ([57]) is one 



component of a system of Roy- type equations [37, 38]. This system was expressed in ref. 



Jl5| ] in terms of the two scattering lengths a^ 2 , a^ 2 . As we do not attempt to solve the 
full system here, we find it more convenient to use <?q(0) as subtraction constant. We 
have included D waves as well into the calculation in order to extend the validity of the 
evaluation of A(t) somewhat above one GeV. The kernels needed in eq. (^) are 

K 00 (s',t) = l (s',t)-I (s',0), 

with I (s',t) = Arcth-^ 

1 {t - Aml){t - Am 2 K ) 2s ~ 2 ^K + t 

Ojf \ Of 

K 10 (s',t) = 3 I (s',t) ( 1 + ^^tt - -^—^ - W,0) 



Kk(s')J Kk{s') 
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K 20 (s',t) = 5 



r f / , ( 6s't 6(s't) 2 
Io(s',t) 1 + - — + 



Kk(s') \l K {s>) 



+ , 2 ; A (s / (-2 S , + 2S^-t) + -(t-4m2)(t-4m^)) -l (s',0) 



G 20 (t',t) 



16 x 5 + i - 4E nK ) 



^3 t / (t / -4m2)(t'-4m^) ' 
The various contributions and the result for A(i) are displayed in Fig. 2. 



(59) 
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Figure 2: Left-hand cut function, A(t), and its various contributions. 



In order to solve equation (^) we first construct the Omnes function over the range 



= exp 



t f*o (f) S (t')dt' 



vr h m \ t'{t> - t) 



= tt R {t) exp[i<f> s (t)e(t - 4ml)6(t - t)] (60) 



where £ln(t) is real. When t approaches the KK threshold, O^(i) has the following 
behaviour, 



lim VL R (t) ~ \t — Am K 



(61) 
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Therefore, two cases must be considered depending whether (j}g (Am 2 ^) is smaller or larger 
than 7T. Let us first consider the case 

M^k) < * ■ (62) 
The solution to eq. ( |57| ) is obtained by introducing the function 

f(t) = ^ r) (g° (t)-A(t)) (63) 

and noting that it is analytic except for a right-hand cut. It can thus be expressed as a 
dispersion relation which is defined up to a polynomial which depends on the behaviour 
at infinity of f(t) |35|] . We will assume that f(t) is bounded at infinity by a polynomial of 
degree one and thus, g® can be expressed in terms of two subtraction constants, 

o/ s s[ „ t 2 f 4m K A( t') sin 6 s (t') 



(64) 



«o = A) = 7^ " A '(°) ( 65 ) 



One observes that the integrals converge at t' = Am 2 K if the condition (|6^) is satisfied. A 
small calculation shows that at t = Am 2 K the condition (g^) output = (9o)input is automati- 
cally satisfied in eq. (^). Concerning the parameters a^, Pq it is not difficult to see that 
for the purpose of using gg with 0(p 4 ) precision it is consistent to use the values of cto, (3q 
with 0(p 2 ) precision, i.e. 

V3_ 

647T/2 

and g® gets fully determined (the value of the derivative A'(0) is determined numerically 
to be A'(0) ~ 0.256 GeV" 2 ). The influence of A(t) is illustrated in Fig.3 which compares 
the full solution from eq. (|M|) to the solution with A set equal to zero. One sees that in 
the energy region where we really need to use eq. (|64]), i.e. below the KK threshold A(i) 
actually has a rather small effect. The solution is essentially controlled from the chiral 
constraints at t = and the experimental input at t > ^m 2 K . Above the KK threshold, the 
agreement of the experimental data with the output from eq. ([54]) seems much improved 
if Ait) is included. In the case where 

5 (4m^) > vr , (66) 

we need only modify the definition of f(t) (eq. (|63[)) to 



f(t) = t - 1 ^(g (t)-A(t)) (67) 

and make the corresponding change in the preceeding formulas. Because of the extra 
factor of t one needs to introduce one more subtraction, and an additional parameter, 70 
appears in the solution: we simply fix 70 so that input-output agreement is retained in 
the physical region when the condition (|6^) holds. 
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Figure 3: Solutions of eq. (57) with A(i) included and A(i) ignored. 



There is a subtlety in the above calculation which must be discussed. Clearly, because 
of the singular behaviour of the term 1/Q#(£') in the integrand at the KK threshold the 
result will be rather sensitive to the value of |<7o(i')l m this re gi° n and, in particular, to 
its value exactly at the threshold, while experimental information starts slightly above 
the KK threshold. A simple solution is to use to slightly larger than 4mj^ in eq. (|6C|). 
We have done so and found good stability of the result. We have also used the following 
simple method which, at the same time, provides an alternative extrapolation method in 
the whole region of interest. We first construct an Omnes function over a region [4m^,ii] 
with ti >> <im 2 K , 

~t_ /■*! Mt'W 

7T J 4ml *'(*' - *) 



fii(£) = exp 



(68) 



and then consider the function 



v °(t) 



Oi(t) 



9°o(t) ■ 



(69) 



The function Vq is analytic with a left-hand cut, and a right-hand cut which only starts at 
t = t\. Therefore, V^ is expected to be a smooth function in the region [0, t\\ and we can 
use approximations by polynomials there. In practice, we used fourth order polynomials, 



V° 
1 o 



(t) = a + Pot + j t 2 + 5 t 3 , 0<t<t fit <h 



(70) 
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The first two parameters are determined from ChPT as above and we fit the remaining 
two to values of V§ determined from the data above the KK threshold. The energy range 
in which the fit is performed t < tfa cannot, of course, be made two large otherwise higher 
order polynomials would be needed. Fig. 4 shows two different fits and illustrates that 
this procedure, while simple, is also quite stable. This procedure allows one to determine 
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Figure 4: Polynomial fits of g® with the right-hand cut (partly) removed (see eq. flgSp). 

the value of |<?o(4m^)| (which is thus correlated with the value of cj)s (4m^-)) and can be 
used for correctly computing the integrals above (64). We can also use the polynomial 
approximation to V® to extrapolate g$ below the KK threshold (note that this method 
requires no knowledge of the left-hand cut and no assumption concerning asymptotic 
behaviour). We found that the two methods of extrapolation are in very good agreement. 

The solution for g§ has a rather strong dependence on the value of the phase (j>s at the 
KK threshold as is shown in fig. 5: the larger </>5(4m|f) the higher is the corresponding 
/o(980) resonance peak. Another source of uncertainty in this calculation is the fact 
that the two available data sets for |<7g|, while having small error bars, are not exactly 
compatible. The data of ref. |28| lies systematically below the data from ref. [^, p0| . 
The corresponding influence in the /q(980) peak is shown in fig. 6. 



Let us now turn to the I = 1 amplitude g\{t). We will again here rely on the exper- 
imental data from ref. p3] above the KK threshold and chiral symmetry at t = 0. Let 



18 



2.5 



£ 1.5 



0.5 



<t> s (t )=150 
* s (to)=175 
<£ s (t )=200 
Etkin 



~i r 



~i r 




Figure 5: Comparison of several solutions of eq. (5j_) for g^{t) corresponding to different 
input values of the threshold phase ^s^m^). 



us first consider the phase of g\ which we denote by (j)p(t). In the range Vt < 0.82 GeV, 



(f)p(t) is equal to the I = 1 irir phase-shift and we use the parametrisation of ref. [ 32 1 
which is constrained from the Roy equations. In the range \/t > 2m k the measured phase 
has been shown in ref. [28] to be well approximated by that of a Breit-Wigner tail of the 
following form 



, _ ,„^G^t)G K {t)/2 
<!] • m*-im p (G*(t) + G K (t)) 1 ' 



with 



m T l + R 2 ql ~ m T 1 + R 2 c 2 



Gir(t) = P 3 P 2 % > G K{t) = P 3 P ' ^2 (72) 



and 

G n = (q 3 jVt)G n , G K = (q 3 K /Vt)G K , q 2 = m 2 p /4 - ml, R = 3.5 GeV" 1 . (73) 

The phase from this formula departs from the measured one above 1.6 GeV but we will 
ignore this discrepancy as in this region the I = 1 amplitude plays little role. Finally, in 
the intermediate region 0.82 GeV <Vi< 2m K we use the interpolation formula 

tan 0p(t) = (a + bt)/(t -m 2 p ), 0.82 GeV < y/i < 2m K . (74) 
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Figure 6: Comparison of two solutions for 
al. 



\n°\ 

\9o\ 



based on the two data sets of Cohen et 
respectively. Above the KK threshold the curves are 
piecewise polynomial fits to the data and below the threshold they are obtained from eq. 



28] and Etkin et al. 



(64). 



From this phase we can construct an Omnes function 



Qp(t) = exp 



t f°° (pp(t')dt' 
7T J U t'{t'-t) 



(75) 



The magnitude of g\ remains to be discussed. As in the case of (/(J we expect that it can be 
expressed with a good approximation as a low order polynomial times the Omnes function 
in the whole energy range of interest. In fact, earlier studies based on extrapolations away 
from the left-hand cut have shown that a constant polynomial is sufficient below one GeV 
[ |39| , EJ. With this in mind, we made a fit to the data between 1 and 1.5 GeV with a 
polynomial containing a constant term plus a term quadratic in t, 



g\(t) ~ «i(l + /3it 2 )ft P (i), Vt < 1.5 GeV 



(76) 



fixing ct\ = v / 2/487r/ 2 from 0(p 2 ) chiral symmetry. A good fit is obtained in this way 
with (3\ = — 0.187GeV -4 such that the quadratic term is indeed small below 1 GeV (in 
discussing the errors we will introduce a linear term as well). The result for \g\ \ is shown 
in fig. 7. 
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Figure 7: Magnitude of the I = 1 partial wave g\ (in units of GeV 2 ) from the construction 



described in the text, compared to the experimental data t 28]. 



We have also included higher partial waves with 1 = 2, 3, 4. For ! = 2 we include 
the resonances /2(1270), /2(1425), /2(1810) with Breit-Wigner functions analogous to eq. 
( pfTj ) and parameters fitted to the data of ref. |3(|. For / = 3 we include p3(1690) and for 
I = 4 the _/4(2050) resonances. In both cases we take the tttt and KK partial widths from 
the PDG flf] 



4.3 Asymptotic region 

Beyond the energy region where the amplitudes are effectively measured in experiments 
we can hardly make better than qualitative estimates. For this purpose we will assume 
that the resonance region matches to a region where Regge behaviour prevails. More 
specifically, we will consider the dual-resonance model for the irK amplitude (see e.g. 

) 



F ± {s,t,u) 
V K * P {s,t) 



-- -\[V K *p{s,t)±V K *p{u,t)} 
T{l-g K *(s))T{l-a p {t)) 
r(l - a K *{s) - a p (t)) 



(77) 
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For the Regge trajectories we take 

a p (t) = 0.475 + ait, ax*(s) 



0.352 + ais , ai = 0.882 GeV" 



(78) 



and for the parameter A we take A = 1.82 which realises an approximate matching to 
the region known from experiment around y/s = 2 GeV. In taking asymptotic limits in 
formula (ff?]) an ie prescription is understood, for instance, s — ► oo means \s\ — > oo and 
s = \s\ exp(ie). One then finds the well known Regge behaviour, 



ImF ± (s,t,u) 



ttA 



*oo, t fixed 



T(a p (t)) 



,a p (t) 



In the case of F + (s,t) one needs to include additionally the Pomeron, 

ImF + (s,0)p mner on = ~CTS , 
7T 



(79) 



(80) 



in which we take a = 2.5 mb (see the discussion in ref. |3^]). We also need ImF ± (s,t,u) 
in the regime where t — > oo and s — » in the integrals G^(n). The model ([FT]) gives a 
Regge behaviour associated with the K* trajectory, 

ttA 

->oo,s fixed ^ 



ImF ± (sA,u) 



(ait) 



a K *{s) 



T{a K *{s)) 

Finally, we need ImF ± (s, t, u) in a regime where s — ► oo and u 



(81) 



in the integrals H b 



n 



In this case, the term VK* P {u,t) in eq. ( jTq) makes no contribution to the imaginary part 
(this, of course, reflects the exact degeneracy of the K* and -fC| trajectories in this model) 
and the term Vx*p(s,t) becomes exponentially suppressed (this term is the amplitude for 
the reaction tt + K~ — ► tt + K~ and the corresponding u— channel is tt + K + — > ir + K + , which 
is exotic) . The influence of these asymptotic contributions can be appreciated from table 
1 below. 





#o + (2) 


#o + (3) 


H^(2) 


^ + (3) 




Gt(3) 


[A 2 ,4GeV^] 


8.05 


4.66 


4.56 


3.11 


1.26 


0.92 


[4GeV^,oo] 


8.75 


0.72 








1.17 


0.13 




Ho -(2) 




H b (2) 


#f(3) 


GT(2) 


GT(3) 


[A 2 ,4GeV 2 ] 


5.81 




2.32 


1.74 


1.76 


1.17 


[4GeV z ,oo] 


2.18 










1.25 


0.14 



Table 1: Results for the high-energy integrals (see eqs. (18) 
of GeV, showing two different integration regions. 



|) ) in appropriate powers 



4.4 Results 

Let us first perform some simple checks. Equating two dispersive representations of 
F~(s,t) we obtained eq. (^): using this relation at t = gives one relation among the 
building blocks of the sum rules 

2i a + H (2) = Ul + G 6 " (2) + ff" (2) (82) 
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We expect some uncertainty because of the relatively slow convergence in the integrals 
Hq{2), GT(2) but there is some amount of cancellation of these effects. Using the phe- 
nomenological input as described above, we obtain, 



2A 1 + Hq(2) ~ 33.17 , mi + G^(2) + F 6 ~(2) ~ 32.07 [ GeV" 



(83) 



Clearly, there is a very reasonable degree of agreement. This provides a check on the 
construction of the nir — > KK P-wave. We have also redone with our input the calculation 



of Karabarbounis and Shaw |M2] which gives the difference in the scattering lengths a, 



1/2 



3/2 



77V (a, 



1/2 




3/2> 



3m, 



8^(777^ + m,K 



■r(m^0) ~ 0.22 



(84) 



to be compared with the result [42] m n (a 



1/2 3/2 



0.26 ± 0.05. We have a comparable 
uncertainty due to a large extent to the asymptotic contributions. A related quantity is 
the polynomial parameter (3~ (eq. fl36|)) for which a rather precise value is predicted by 
the chiral expansion, 

Sir = 0.25 + 0.01 + 0(p & ) , (85) 

where the successive contributions are shown. Using the sum rule expression (54) for j3~ 
we obtain 



fl(3- 



0.24 



(86) 



which is within 10% of the result from ChPT. The size of the uncertainty in this calculation 
is approximately 15%, so the agreement is satisfactory but it is not possible to separate 
the purely 0(p 4 ) part (in other term, we have no sum rule for L5). We note also that the 
0(p 4 ) contribution being suppressed, the 0(p 6 ) one could be of comparable size. 

Let us now discuss the results for the chiral couplings L±, L2, L3 



these are obtained by first generating sum rules for the polynomial coefficients , /\ 2 , 



We recall that 
^2 

X 1 : here the contributions from the asymptotic regions are suppressed so we can expect 
rather good accuracy. Our results are collected in the last line of table 2, they complete 
and update those already givenQ in ref. [15]. The way in which the errors quoted in the 



10 3 x 


Li 


L 2 


L 3 


L 3 + 2Li 


7T7T ST 0(p i ) 




1.02 ±0.05 




-2.78 ±0.32 


K u 0(p 4 ) 


0.46 ± 0.23 


1.49 ±0.23 


-3.18 ±0.85 


-2.26 ±0.97 


ttK sr 0(p 4 ) 


0.84 ±0.15 


1.36 ±0.13 


-3.65 ±0.45 


-1.97 ±0.34 



Table 2: Results from the sum rules for the chiral couplings L\, L2, L% (multiplied by 10 3 j 
at the scale \x = 0.770 GeV (last line), compared to the results from nir sum rules and 
from the form-factors. 

table are evaluated will be explained in more details below, they do not include any effect 

1 In ref. a factor f%fj< was used in eqs. ( fl7| ) (S) instead of Here, we prefer not to include 
incomplete parts of the 0(p e ') contributions. 
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from 0(p e ) corrections. The table also shows for comparison the results obtained from tttt 
sum rules (taken from ref. [34] in which 0(p 4 ) matching is used) and the results based on 



the K14 form- factors also computed at chiral order p 4 . The numbers quoted in the table 
are taken from the fit of Amoros et al. |43j based on the data f\ of Rosselet et al. [44]. 



We note that the value of L\ obtained previously in ref. ||, 10 3 Li = 0.65 ± 0.27, while 
compatible, has a somewhat larger central value. Our results for L2 and L3 agree within 
approximately 30% with the results from K14 or tttt. Concerning L2, the difference beween 
the K14 and the tttt result is substantially larger than that, our result happens to lie in 
between these two. A discrepancy at the 30% level can be expected as a consequence of 
unaccounted for 0(p 6 ) effects. In the case of L%, however, there is a larger discrepancy, by 
about a factor of two, between our value and that from ref. [43]. The combination L2 — 2L\ 



is suppressed in the large N c limit || . We indeed find a suppression of the value of L2 — 2L\ 
(compared, say, with L\{m p ) or L2(rn p )) even though this combination is dominated by 
scalar resonances. A calculation of the from factors in ChPT to order p e was recently 
performed ]^3|, and the couplings L\, L2, L% were then redetermined, using a model 
to estimate the 0(p 6 ) couplings Cj(yu). The following numbers are obtained |43[ 

10 3 Li = 0.53 ± 0.25, 10 3 L 2 = 0.71 ± 0.27, 10 3 L 3 = -2.72 ± 1.12 [K H , 0(p 6 )] . (87) 

Clearly, variations larger than naively expected can occur as compared to the 0(p 4 ) deter- 
mination. It remains to be seen, and this would be an interesting check of the convergence 
of the SU (3) chiral expansion, how the differences in the results from the various methods 
of determining L\, L2, L3 are reduced, once the 0(p e ) contributions are included. 

In order to estimate the errors, firstly, we have varied the S and P wave ttK phase- 

1/2 3/2 1/2 

shifts inside bands of half-width A5 = 2°, A5 = 1°, A<5 1 = 1°, which correspond to 

the average experimental errors in the region of elastic scattering (which makes the most 

important contribution). For the tttt — ► KK P-wave, g\, we have varied the coefficients 

of the normalising polynomial (see eq. (|76|)), allowing for a term linear in t, and requiring 

that the \ 2 does not exceed twice its minimal value. The coefficient a\ is kept fixed since 

its variation can be considered as an 0(p 6 ) effect, which we do not try to estimate. This 

procedure generates a variation of the height of the p resonance peak in the 10% range, 

which may seem rather small, but affects the results quite substantially, as can be seen 

from table 3. For the S'-wave, g$, we have varied the phase at threshold <j>s(km 2 K ) (which, 

as we have seen, is the parameter on which the size of the /o(980) peak mostly depends) 

in a range between 150 and 200 degrees. Above threshold we have made \g$ \ to vary in the 

whole range allowed by the two incompatible experiments. We have also allowed a 10% 

variation of the scattering lengths a]J 2 a^ 2 (keeping, however, the difference fixed) and, 

finally, in the Regge region we have assumed a 100% uncertainty. We show the individual 

impact of these variations on the sum rule results in table 3. 

We now come to the discussion of L4. From relation (^0|) an important remark can 
be made concerning the convergence: while j3 + and (3~ , separately, contain integrals 

2 We thank P. Talavera for communicating the values of the errors corresponding to this fit. 

3 An indicative fit using preliminary data from the more recent E865 experiment p3(] is performed in ref. 
p^ ] which gives essentially the same central values for Li, L2, L3 and error bars reduced by approximately 
a factor of two. 
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0. 
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3.6 




0.1 


1.5 


2.5 


0.6 


1/2 

a 


0.3 


0.4 


1.4 


0.6 


9\ 


9. 


9. 


36. 


10. 




4.5 


0. 


0. 


15. 


Regge 


0.3 


1.6 


1.0 


0.1 



Table 3: List of different sources of errors (see text for details) and their impact on the 
determination of the Li 's. 

which are slowly convergent at infinity, L4 involves the difference, which has much better 
convergence properties. Indeed, consider the high-energy contribution 

IP + -P-]he = G+(2)-G^(2)-(H+(2)+H~(2))+2^ K (H + (3)-2H+(3)-H-(3)). (88) 

On rather general grounds, the leading Regge contribution in G^(2) and Gr(2) is the 
same and will cancel out in the difference. Also in the second potentially dangerous term 
H^{2) + H^(2) the relevant cross channel is pure 1 = 3/2 and has no leading Regge 
contributions. The other terms in eq. ( |88| ) are, as we have seen, rather insensitive to the 
asymptotic region. Therefore, we expect the uncertainty in L4 coming from the asymptotic 
region to be small, of the same size as in L\, L2, L3. Let us now consider the energy region 
below 1 GeV. Using eqs. (|5C|) and (|54"|) we can write, 



[U]le = 2Af + 2A\ - B+ - Bo +v - Ul ,+2E^ K (-C+ - C ~ + w ) (89) 

(where we have used = A^ , A\ = A\ which is true up to 0(p 8 )). This expression 
contains P-wave contributions, which may seem surprising in view of the well known 
resonance saturated expression M 



jres _ c dC m CdC m 

--wJ+mT, (90> 



which involves only scalar resonances. It is in fact possible to write an alternative expres- 
sion for (|89| ) : using the crossing symmetry relation (|8^) the P-wave combination 2Af — u\ 
gets replaced by contributions from above the resonance region and the P-wave term Af 
has, in fact, no contribution from the resonances. This alternative expression has only 
S-wave resonance contributions but is not as rapidly convergent. 

Numerical results for L4 are shown in table 4 for several input values of the threshold 
phase </>s(4m| i -) . We have also mentioned that the two experimental measurements of 
|(7q| of refs. [p8| and [ 30 1 have a somewhat inconsistent normalisation. We have performed 



the calculation for each data set separately. A clear feature from these calculations is 
that L\[rrtp) is suppressed, and has a magnitude similar to 2L\ — Li ■ It is not very 



25 





150° 


175° 


185° 


200° 
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0.08 
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0.18 
0.10 


0.22 
0.13 


0.27 
0.16 


0.34 
0.21 


Etkin 
Cohen 



Table 4: Sum rule results for L^(p, = 0.770). 



easy to decide on which central value to choose. We will make the choice of believing 
the data of Cohen et al. |28| for the value of (/>s(4m^-) , which is then close to 200°, as 
they argue that the presence of a P-wave in their experiment (which is absent in the other 
experiment) helps in correctly determining the S-wave phase at threshold. Concerning the 
normalisation of \g$\ above the threshold we may average over both experiments. Taking 
into account the main sources of uncertainty (see table 3) we would then obtain, 

L 4 (a« = 0.770) = (0.22 ± 0.30) • 10~ 3 . (91) 

This result is a refinement of the previous estimate of Gasser and Leutwyler Q, L^fj, = 
0.770) = (-0.3 ± 0.5) • 10~ 3 , based on the assumption that OZI suppression holds but 
without precisely knowing the value of the scale at which it does. Other results can be 
found in the literature pq, 47] which, however, are based on some assumptions allowing one 



to determine scalar form- factors. The fact that they agree with (|91| ) indicates that these 
assumptions are reasonable. Finally, Amoros et al. [p| have attempted to determine 
L4 from Ki± data, using their 0(p 6 ) calculations, and they obtain L 4 (/x = 0.770) = 
(—0.2 ± 0.9) • 10~ 3 which, as expected, is not very tightly constrained. 

5. Conclusions 

In this paper we have performed an evaluation of the set of sum rules proposed in ref. 



[pup. While some results were already presented in ref. IS], the calculations performed 
here are more complete: they take into account the contributions from partial waves 
beyond the S and P waves as well as asymptotic contributions. In order to fully exploit 
these sum rules, one needs to perform an extrapolation of the I = and I = 1 partial 
waves of the tttt — * KK amplitude. This can be performed using standard Muskhelishvili- 



Omnes techniques and was considered a long time ago [16]. A great improvement over 
these calculations is the availability nowadays of direct and precise experimental results 
concerning the tttt — ► KK amplitude, so there is no need to make use of inelasticity in tttt 
scattering, which is not very precisely determined, and requires the assumption of exact 
two-channel unitarity. We have checked the stability of the calculation by comparing 
different approaches to the solution. The main source of uncertainty comes from the value 
of the phase at the KK threshold because the height of the /o(980) peak is strongly 
correlated with the value of this phase. The two available experimental data sets do not 
agree on this value and we have assumed a plausible range of variation. We have obtained 
a redetermination of the three 0(p 4 ) chiral couplings L\, L2, L%. Comparison with former 
determinations allows a test of the SU(3) chiral expansion. For instance, it is encouraging 
that the value of L2 that we obtain is intermediate between the value from K^ and the 
value from tttt sum rules. Besides, we have obtained for the first time an evaluation of L4 
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at the same level of precision and reliability as L\, L2, L3. This was not possible from the 
K14 from-factors because of an accidental suppression of the coefficient of L4 in this case. 



In constructing the I = and I = 1 partial waves of the irir — > KK amplitude, we have 
solved a subset of the system of Roy-Steiner equations. A further improvement, which we 
have not performed here, would be to use the full set of equations in order to constrain the 
low energy part of the ttK — > ttK amplitudes. We note however, that the range of energies 
where this is needed, that is, between the threshold and the energy where the data start is 
smaller for ttK than it is for -kit. Solving these equations would help in deciding whether 
a strange counterpart of the a meson, the n meson, actually exists (e.g. ref. |0] and 



references therein) or not [48|. An obvious further improvement would be to use the sum 
rules in association with a chiral 0(p e ) calculation of the irK amplitudes. These, taken 
together with the mr sum rules (associated also with an 0(p 6 ) SU(3) calculation of the irir 
amplitude) and the available calculation of the Kn form- factors to this order [45 would 
no doubt greatly improve our understanding of the chiral expansion in m s . 
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